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We revisit the expression for the conductance of a general nanostructure - such as a quantum 
point contact - as obtained from the hnear response theory. We show that the conductance 
represents the strength of the Drude singularity in the conductivity a{k,k';iuj 0). Using 
the equation of continuity for electric charge we obtain a formula for conductance in terms of 
polarization of the system. This identification can be used for direct calculation of the conductance 
for systems of interest even at the ab-initio level. In particular, we show that one can evaluate 
the conductance from calculations for a finite system without the need for special "transport" 
boundary conditions. 
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I. INTRODUCTION 

The conductance gives the current as a linear response 
to an applied voltage for a given finite sample, specified 
by its atomic geometry. Typically, the sample of interest 
may be a thin chain of metallic atoms or a single molecule 
attached to metallic electrodes. The usual theory for 
ab-initio prediction of conductances for these nanoscale 
devices is based on the Landauer formula^ 



conductance formula^ 



2e2 
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which identifies the conductance with the transmission 
probability for an electron at the Fermi level, T^ep), to 
penetrate the sample. 

Given the impressive success of the Landauer-Biittiker 
equations'^ in many mesoscopic phenomena, many au- 
thors tried to relate it to basic principles. Since we are 
concerned with linear response, the general Kubo lin- 
ear response formalism^ is well suited to this purpose. 
Among the first to do so were Economou and Souk- 
oulis^ who presented a derivation for noninteracting elec- 
trons. This was subsequently vigorously discussed^iSii, 
since perfect transmission T = 1 results in a finite con- 
ductance, not infinite as might be expected. Eventu- 
ally the debate was reconciled with the realization that 
that Formula Q gives the conductance between reser- 
voirs and =|- represents an unavoidable minimum con- 
tact resistance if such macroscopic electrodes are con- 
nected through a quasi-lD contact. However, there are 
no such electrodes in the Kubo-like derivations, as has 
been pointed out by Landauer^, and so the result re- 
mained somewhat puzzling. The argument in favour of a 
Kubo-like derivation was an introduction of an auxiliary 
external field acting in the leads that compensates ex- 
actly the charge-density oscillations produced by the ho- 
mogeneous field in the sample, resulting in the "4-point" 
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that relates the current to the local drop in total elec- 
trostatic potential. However, in many experiments and 
situations it is the 2-point conductance G that is re- 
quired, in which case Formula is frequently used for 
conductance calculations (e.g. for nanowires and metal- 
molecule- metal junctions). Therefore, clarification of the 
relationship of this formula to the Kubo approach re- 
mains necessary, and may lead to practical advances for 
ab-initio calculations. 

Much of the confusion arising from the "boundary and 
bulk contribution to the conductivity" introduced by the 
ideal leads has been clearly resolved by SolsS who em- 
phasized the importance of global charge conservation 
and its relation to gauge transformations. In that paper, 
however, only the result ^ was derived and discussed 
and the induced charge-density oscillations were ignored. 

The fundamental understanding of conductance as 
quantum-mechanical transmission relies on the picture 
of distinct electro-chemical potentials for the two elec- 
trodes. These dictate the occupation of states in the 
non-equilibrium sample. An alternative point of view 
(equally valid but much closer to the Kubo formalism) is 
that of electrons accelerated in the applied external field. 
The ability of the field to accelerate electrons in extended 
or bulk materials has been used by KohniS to distinguish 
between conductors and insulators where conductors are 
characterized by \im^^oQ{a{ui,q = 0)} ~ A/u},A ^ 0, 
and the case A = characterizes insulators. This point 
of view, but applied to conductance of finite samples or 
mesoscopic conductors, has been considered by Fenton^i 
who strictly differentiated between a localized external 
field leading - through acceleration of electrons over a 
finite distance - to a finite conductance and yielding the 
Landauer formula, and a homogeneous field giving a re- 
sponse characteristic of the bulk material with lu~^ sin- 
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gularity. An approach based on coupled self-consistent 
equations for the current density, electron density and in- 
duced field, implicitly also representing the "conductance 
based on acceleration" point of view, was set out in de- 
tail by Kamenev and KohniS. They focus on the 4-point 
conductance, which they define as the current divided 
by induced potential drop, thereby completely neglecting 
the electrostatic drop in the external potential. This is 
appropriate in some cases, since a finite constant exter- 
nal field over infinitely long leads contributes negligibly 
to the total potential drop around the sample. It is im- 
portant to note that their self-consistent field is obtained 
from the ID Poisson equation, i.e. a system resembling 
parallel planes of charge e moving along their normal 
rather then electrons confined to a ID wire. The former 
systems has far stronger screening behaviour then the lat- 
ter. More detailed calculations by Sabliko\>i^, who also 
used the Kubo formalism, show that the effective interac- 
tion between electrons in a ID wire is not strong enough 
to screen the charge and a uniform current coexists with 
an induced dipole moment of the electron density in the 
wire. 

In this paper we derive a formula for the steady-state 
conductance in terms of the irreducible polarization of 
the system of interest. Our starting point is application 
of an external field and following the time evolution for 
large times. We will show that only the q = Fourier 
component of the field influences the steady-state current 
because the nonlocal conductivity is a sharply peaked 
function of q, eventually becoming a Dirac delta func- 
tion. The weight of this delta function has directly the 
meaning of conductance. At the same time, this singu- 
lar behaviour also corresponds to the Drude singularity 
3{(t(u;)} ^ uj^^ if we set q = and then take the limit 
Lu 0. It is interesting to note that from this point 
of view we can identify the conductance with the Drude 
weight, a connection encountered in somewhat different 
circumstances in the conductive properties of extended 
systemsiii. In addition to providing an interpretation 
of the conductance, our formula is suitable for practi- 
cal ab-initio calculations. Using the continuity equation 
we circumvent the need for matrix elements of current 
operator, which are in general difficult to evaluate nu- 
merically, and express the conductance in terms of the 
irreducible polarization at small imaginary frequency - a 
quantity that can be evaluated efficiently even using ab- 
initio methods that describe elecronic correlationi^ii&. 



II. CONDUCTANCE IN IMAGINARY 
FREQUENCY 

Consider a system, infinitely long along the x direc- 
tion and finite, infinite or periodic along the other two 
perpendicular directions. Along the x axis we apply an 
infinitesimally weak external electric field. The response 
of the system will be in general non-local in time and 
space, and quantitatively described by the Kubo linear 



response theorji^: 

j(r,i) = f a{r,r'-t^t')-^{r',t')d\'dt', (3) 

J — OO 

l—ill 

a{r,r';t) = ie{t) / drTr {p,,j(r', t + r)j(r, 0)} (4) 
Jo 

where j(r,i) is the current density, E{r,t) external elec- 
tric field, CT(r, r';t) non-local tensor of conductivity, 
Q{t) = I for t > and Q{t) = otherwise and p^q 
the equilibrium density matrix. 

We take the form of the field as Ex{r,t) — 



E{x)e{t),Ey 



0. Introducing the current 



I{x,t) = J dS-i{r,t) (where we integrate over the surface 
perpendicular to x) and utilising our particular choice of 
external field we obtain 



I{x,t) = J dx'I„{x,x';t)E^{x') (5) 
h{x,x';t) = I dt' f dS -air, r';t')-dS' (6) 

J —OO J 

In the steady state we are interested in the limit t — s- oo 
so that for the steady-state current we have 



Ic,{x,x') 



lim / dt'e-°'*'a{x,x';t') 

<^^"+ J-oc 

a{x, x'; Lu — iO^), 



(7) 



where we have introduced the effective one-dimensional 
conductivity cr(a;, 2;'; w) = / dS • CT(r, r'; oj) • dS'. The pos- 
itive infinitesimal a plays an important role for systems 
without dissipation, where application of a field for an 
infinitely long time results in the system being heated. 
This is avoided by introducing an effective finite (but 
large) measurement time T ~ a~^. In this way we avoid 
the heating problem because we first let E{x) — > (linear 
response) and only afterwards a 0~^. (Identical ex- 
pressions are obtained assuming adiabatic switching-on 
of the external field E{x, t) — E{x)e°'* as t approaches 
from below, and measuring the current at t = 0~, which 
is the more conventional point of view.) 

We are principally concerned with the total current / 
and not the current density. Similarly, we would like to 
work with the bias voltage V ^ j dr ■ E(r) — J dxE{x) 
rather than the field itself. To achieve this we Fourier- 
transform with respect to x 



dq' 
2^ 



a{q,q')E{q') 



a{q,q') = / e~*«^+*«'=^'a(a;,a;')da:dx' 



(8) 
(9) 



First, the steady-state current needs to be independent 
of X, I{x) = /, as a consequence of the equation of con- 
tinuity, so that I{q) = 2TrId{q) and therefore directly 
also a{q,q') ^ S{q). Second, as a direct consequence^ 
of the linearity of the theory, the steady-state current is 
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uniquely given by the bias V = E{q = 0) only. This 
means that for two external fields E{q) and E'{q) whose 
long-range parts are equal {E{q = 0) = E'{q = 0)), 
but otherwise are arbitrary, one has to obtain identi- 
cal steady-state currents. We therefore conclude that 
c('Zj'?') ~ ^W) ill the limit a — > 0^, i.e. it is a sharply 
peaked function of q' so that, based on Eq. (jHJ, it is only 
the values of E{q) or E'{q) at q = that affect the steady 
current. Using these observations in Eq. ((SJ we obtain 



dqdq' 
47r2 



(10) 



or using the fact that a{q,q') 
write 

lim a(q,q';ia) = 4Tr'^G5(q)5(q') 

Q^0 + 



5{q)6{q') we can also 



(11) 



This is one of the central results of this paper and directly 
shows that conductance is the strength or weight of the 
Drude singularity. Formula (|10|l can be used for extrap- 
olative evaluation of the conductance from the conductiv- 
ity at small imaginary frequencies, and since it is peaked 
function in q, g', the conductivity needs to be evaluated 
for small values of these variables only. We will demon- 
strate this approach numerically in Section Ivl 



III. CONDUCTIVITY AND POLARIZATION 

Evaluation of the conductance using Equation 1)10(1 re- 
quires calculation of the Fourier transform of the con- 
ductivity. The latter is, in a real-space representation, 
given by Equation @ which itself is difficult to evaluate 
for all necessary r,r'. On the other hand, substantial 
experience has been accumulated in calculations of the 
polarization x(r, r';t) , defined hyi^^ 



X(r,r';i-t') 



Sn{r, t) 
SV{r',t')' 



(12) 



where Sn{r,t) is a change in density due to infinitesi- 
mally weak external potential V{r,t). To relate these 
two we utilize the equation of continuity integrated over 
the cross-sectional area 



d^Iix,t) = dtN{x,t). 



(13) 



From the definition of the polarization we have N{x, t) = 
dt' J dx'x{x,x';t — t')V{x,t'), where our x(a;,a:';t) is 
now, similarly to a{x,x';t), integrated across the cross- 
sectional area (dS and dS' integrals in ((HJ). Substituting 
the linear response formulae for current |SJ| and density 
N{x,t), Fourier-transforming into q,q' variables and us- 
ing the fact that the external potential is arbitrary, we 
immediately obtain 



a{q, q'] ia iO+) = lim — x(g, ia). 

a^a+ qq 



(14) 



We note that the singular character of a{q, q') for small 
q, q' as given in Hill) does not arise fundamentally from 
the 1/qq' prefactor, since x(9, '^) ~ Q<l'fil,Q'','^) where 
/(O, 0) 7^ 0. This property is a consequence of conserva- 
tion of total number of particles {n{q = 0) = TV) or the 
absence of response in density if we change the potential 
everywhere by a constant q' = 0. 

Expressing the conductance through the polarization 
is particularly suited to a correct treatment of a system 
of interacting electrons. In the first place, it is crucial 
to define the conductance not as a coefficient for current 
dependence on the external but on the total electric field 
E*(r) = E(r) -h E*(r). The induced field E'(r) can be 
obtained from Poisson equation 



iiq± + G) • E'(g, G; t) = ^AirSniq, G; t), 



(15) 



where ic is a unit vector in the x direction and G is a 
2D reciprocal lattice vector corresponding to the perpen- 
dicular coordinates. Using the linear response result for 
Sn{q, G; t) in terms of x(r, t) (see Eq. (|12|l ') we arrive at 
the relation between the total and external field 



E*(g,G) - J dq'[ScSG,oS{q-q' 
qx + G 



(16) 



Att 



q' q^ + \G\- 



x{q,G;q\0) 



E{q') 



If we multiply the latter equation by the unit vector x and 
subsequently invert it for G = 0, we obtain the relation 
between the external field E{q) and the total field along 
X averaged over the cross-sectional area, £'*(q, G — 0). 
Substituting this relation into equation Q we obtain, 
after some algebra (Appendix A), 

G^l^-f ^^M^dqdq', (17) 
a^o+ 47r"^ J qq' 

where we introduce the transport part of the polariza- 
tiori\ X*('7i ^'^)^ related to the irreducible polarization 
)(^{q, G; q', G') through an equation of Dyson type 



X*(g,G;(?',G') =x°(g,G;<7',G' 



(18) 



+ [dkY, X*('Z,G;fc,K)— i^x°(fc,K;g',G'). 

The corrections entering through K terms in equa- 
tion H18() are known as local-field effects in the context of 
evaluation of the macroscopic dielectric functioni^. Here, 
however, it is not quite the same since fc = 0, K 7^ part 
is included into x* whereas for the macroscopic dielec- 
tric function the sum is restricted to fc 7^ 0,K 7^ 0. The 
omission of the K = term involving |f in equation 
(|18|1 stems directly from the fact that it is precisely this 
term that converts the drop in external field into drop in 
the total field, as can be seen from Eq. H16() when tak- 
ing into account that the drop in total field is given by 
Ay* = E\q = 0, G = 0). Correct evaluation of the 
conductance in 3D therefore requires inclusion of these 



4 



"perpendicular" local-field effects, included within x* but 
not in x*^, into account. Essentially, x' describes the re- 
sponse of the density to the effective potential, except 
that long-range screening of the potential in the x direc- 
tion is specifically excluded, allowing the conductance to 
address the applied voltage rather than the local poten- 
tial drop. 

For systems translationally invariant along the perpen- 
dicular directions, x^il-^ 0; k, G) ~ 5o,g and the last term 
in (|18|l becomes identically zero. It follows that in this 
particular case x' = X° and the conductance of non- 
interacting electrons, defined with respect to drop in the 
external field, is identical to the conductance of inter- 
acting electrons, treated within the RPA approximation, 
defined with respect to drop in the total field. 




FIG. 1: Semi-classical interpretation of a finite-field potential 
V{x) appfied to a wire. The local Fermi spheres are depicted 
as for 2D gas for clarity. The filled area on top of the equilib- 
rium spheres indicates accelerated electrons, while the empty 
area inside the equilibrium spheres represents states emptied 
by the deccelerating effect of the uphill field. 



IV. LANDAUER FORMULA 

In this section we present simple, analytically tractable 
cases that illustrate the theory of the preceding sections. 
Consider first the case of non-interacting electrons in a 
quantum wire with only one subband. It is well known 
that the polarization function has the formii 



1 



7, — log 



^2 - (g2/2 - kpqf 



x27rS{q-q'), 



w2 - {q^/2 + kFq)^ 

(19) 

where the factor 2T:d{q—q') arises trivially from the trans- 
lational invariance of the system along the x axis. Using 
this expression in H14|l immediately gives a{q,q';ia 
iO+) = 2TTS{q')S{q) and therefore through Eq. ^ G = 
l/27r, the quantum of conductance (i.e. e^/h in S.I. 
units). Using Eq. H17|l without application of the limit, 
and Eq. (|19|l . it is also possible to obtain the analytical 
dependence of conductance G on imaginary frequency 
a — ^iLu for this system. 



G{a) = 



1 



(20) 



This functional form will be extremely useful for numer- 
ical calculation of G{0~^) based on extrapolation to zero 
frequency, described in Section Ivl 

It is also instructive to explore the behaviour of 1^ {t 
oo) directly (see Eq. ©). The latter can be calculated, 
e.g. using Fourier transform of (|19|) 



'dt'aiq,q';t')^ fdt'^J^^^ 
Jo Ql 

^^x2nSiq-q% (21) 



where we have used the expression (I14f) in the real-time 
domain. The result is 



2 0^ 

sm{kfqt) sin(-^i)e(<) x 2TT5{q - q'), 



nqH 



(22) 



which for t — s- oo evidently approaches the form Hll|) 
with G = The lesson from these analytical examples 
is that the well-known formulae for polarization can be 
used for evaluation of the conductance of the system, 
and that the limit iuj 0+ clearly corresponds to the 
steady-state limit t oo. 

Because the conductivity is sharply peaked around 
q' = 0, only the q — component of the applied field 
E has an effect on the total current. This allows us to 
include, as a limiting homogeneous field. For a 

system with a constant field applied over length L, the 
field itself must be '--^ 1/i so that we have a sufficiently 
small finite drop V = Jj^ Edx ~ 1 . The physical mean- 
ing of the finite conductance of a scattering-free segment 
of a metallic wire is clearly a manifestation of the free 
acceleration of electrons over the distance L. This is in 
agreement with the point of view advocated by Fentonii. 
However, there is no need to keep L strictly finite; in 
fact the limiting case L — > cx) can be still characterized 
by a finite overall conductance by taking careful limits 
and using properties of the Dirac delta function, as men- 
tioned above for the case of the free electron gas (follow- 
ing Eq. (|22|) '). In simple terms, this represents the limit 
of increasing the length over which the external field is 
applied, L, while decreasing the intensity of the field E 
in such a way as to keep the drop V ^ EL finite and 
small. 

A simple qualitative demonstration of this argument 
can be made semi-classically. Consider a finite drop over 
finite region of a ID wire as shown in Fig.^ The current 
can be obtained from the local Fermi distribution on the 
far right. Electrons that are occupying states above the 
equilibrium distribution have travelled to the right from 
the region with an accelerating field. These, during their 
fiight through the region, gained energy ISV so that they 
represent current / = kp where we sum over states 

with accelerated electrons, i.e. k : ek & {eF,eF + AV). 
This obviously leads to the quantum of conductance 
G — While formally we have striking similarity with 
the usual "two chemical potentials" picture, it should be 
noted that the microscopic interpretation of these expres- 
sions is different. As can be seen from Fig. Q] the local 
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charge neutrality is clearly violated in the constant po- 
tential regions on the right and left while in the center 
it is locally charge neutral where the distribution corre- 
sponds to the shifted Fermi sphere. The infinite L limit, 
discussed in previous paragraph, resolves this problem by 
sending the regions with unbalanced local charge away 
into infinities. 

When we insert a localized scattering potential char- 
acterized by the transmission matrix 



T(fc) 



r(k) t(k) 
t{k) r{k) 



(23) 



into fD gas of electrons, the Landauer formula is ob- 
tained. We obtain the desired demonstration of the 
sharply-peaked character of a{q,q';ia) by analytically 
continuing Kamenev and Kohn's expressions^ onto the 
imaginary frequency axis, 



a{q,q';ia) 



1 

X 



2kpa 



kpq^ - 
2kpa 



;2nS{q-q')-\r{ep)\' 



2kpa 



2 „2 



27:\t{eF)\'6{q)6{q'). 



(24) 
(25) 



Alternatively, we can arrive at the same result using 
Equation (|17|) , which explicitly shows the utility of refor- 
mulation of transport through the polarization function 
in imaginary frequency. A detailed derivation is given in 
Appendix B. 

We stress that this result for the two-point conductance 
of non-interacting electrons also applies to interact- 
ing electrons if the conductance is defined with respect 
to the total field at the RPA level of approximation. The 
4-point conductance |2Jl may be formulated in a similar 
fashion, which, if combined with the approximate effec- 
tive electron-electron interaction vl^{q) ~ -i^, would be 

equivalent to the treatment of Kamenev and Kohn^^.(In 
their calculation the 4-point conductance is defined as 
G^^ = I/V^, where is the self-consistent induced drop 
only. The reason for neglecting the drop in external field 
is due to the above mentioned limiting procedure; since 
E X L is finite while i — > oo, i? x L' will be negligible for 
any finite L'. Since the 4-point measurement is meant 
across a finite distance L', the contribution to the to- 
tal drop on finite distance comes solely from the induced 
field.) 



V. FEASIBILITY FOR NUMERICAL 
CALCULATIONS 

The expression for the conductance H17|l , together with 
some experience gained from the evaluation of G for 
lead-scattering-lead system in the Appendix, motivates 
the following suggestion: can the correct conductance 
of an infinite system be calculated using a finite system 
with some arbitrary boundary conditions, determining 



the polarization as a function of a, and then extrapolat- 
ing a ^ 0? 

For the purpose of this exploration we have considered 
a square-barrier potential in ID, for which the results 
for transmission coefficients are well known analytically. 
We calculate the polarization function using the Green's 
function of the equilibrium system given by equation 



H{x) Q{x,x';e) ^ d{x — x'), 



(26) 



supplemented by the chosen boundary conditions. For 
simplicity and for the purpose of illustration, we take 
the extreme case of "zero" boundary conditions on the 
wavefunctionsifi-. In this case, it is clear that G{a) 
as a — !■ 0. In terms of the Green's function we can easily 
express the polarization as 

occ 

x{x, x'; ia) = ^ {g{x, x'; -I- ia)(l)^{x')(l)*{x) 

i 

+ g{x', X- e, - ia)(^^{x)(j)*{x')) , (27) 

where we sum over occupied states only, and 0i are eigen- 
states of H from H26I) . The Green's function at given en- 
ergy is found by direct integration of (|26|) and eigenstates 
are easily found by matching plane waves with the chosen 
boundary conditions. Finally we use a discrete Fourier 
transform, x(a;,a;') — > x{K,K'), to obtain the estimate 
of conductance 



^/ N X{.K, K';ia) 

G{a) — > a . 

^ ' ^ KK' 

K,K' 



(28) 



The discrete Fourier transform introduces a convergence 
parameter: the real-space step size Ax or, equivalently, 
the number of discretization points N. Altogether, the 
calculation needs to be converged with respect to the 
system size L and with respect to N, and extrapolated 
with respect to imaginary energy a. 

Before we discuss results for non-zero barrier height, let 
us first consider the case of zero barrier height, i.e. the 
free electron gas. In this case we know exactly the whole 
dependence of G{a), given by Eq. (|^. We have fixed 
the Fermi wave- vector to kp — 7r/3 a.u., corresponding 
to ep ~ 0.5 a.u. In figure |21 we show numerical results 
obtained using Eq. The labels ^^L-N" represent 

system of length L a.u. with N discretization points in 
real space, or equivalently, the number of K points in the 
discrete Fourier transform. We obtain our chosen Fermi 
energy for lengths L = 48 and 96 when occupying 16 and 
32 states respectively. 

From the figure ^ we see that the numerical results 
converge to the analytical expression for energies a S> 
ep. The lower limit of this range can be brought closer 
to by increasing the length of the system. However, 
for a longer system the convergence with respect to N 
becomes more demanding, with the error growing with 
energy. This behaviour can be traced to the cusp in the 
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FIG. 2: The dependence of conductance of free electron ID 
gas on imaginary frequency. Black line corresponds to the 
analytical result 1201 and the numerical data labels "L-N" 
represent system of length L a.u. with A'' discretization points 
in real space. Apart from the low-energy drop to zero, which is 
characteristic of a finite system, the functions converge to the 
exact result for the infinite system. This opens the possibility 
of calculating the D.C. conductance by extrapolation from 
moderate values of a. 



polarization xi^i Sii x — x' . It can be easily found 
from Eq. 127|l that the size of this cusp is 

occ 
i 

This expression could be used to remove the cusp from 
numerical Fourier transforms and substantially decrease 
the needed number of K points. In the present paper, 
since our system is computationally undemanding, we 
have instead simply increased N until satisfactory con- 
vergence was achieved. 

Next we include a square barrier into our system. In 
Fig. 121 we show results of calculation with barrier with 
width a = 2 a.u. and height ep. In this case the transmis- 
sion probability is approximately 50% which is indicated 
in the upper graph with arrows. In the lower graph we 
show the dependence of conductivity on the imaginary 
frequency a. The two sets of numerical data correspond 
to odd and even numbers of occupied states, converged 
with respect to N within the considered range of a. The 
convergence of data with respect to L is first fast but 
eventually becomes rather slow. In fact, it is not our aim 
to get the conductance of an "almost infinite" system, 
but rather to be able to extract the conductance of an 
infinite system from a calculation for a sufficiently large 
but finite system. The shown length L = 96 a.u. is to be 
compared with the Fermi wavelength Af ~ 6. 

A striking difference between odd and even number of 
occupied states comes from the fact that placing the bar- 
rier of length a = 2 in the center affects the symmetric 
and antisymmetric states (at the Fermi energy) differ- 
ently. Clearly, the former have a much larger amplitude 
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FIG. 3: The dependence of conductance, given by Landauer 
formula, on the Fermi energy {kp = 7r/3 in our calculations) 
for a square barrier of height ep and width A_f/3 (upper 
graph) and the dependence of the numerically calculated con- 
ductance on the imaginary energy (lower graph). G{a) ap- 
proaches the exact result of infinite system before it turns 
rapidly towards zero. 



at the scattering potential than the latter (which are zero 
in the center of the system) and therefore the symmet- 
ric states will be more affected. This is in agreement 
with the conductance for 32 occupied states being below 
that for 31 occupied states. For L ^ oo these two states 
become degenerate and the difference between odd and 
even cease to exist. By bracketing the exact conductance 
in this way, we can obtain acceptable convergence even 
with the crude "zero" boundary conditions. 

When calculating conductance we need to be able to 
extrapolate our data to zero energy. As we have pointed 
out above, the small-energy turnaround of G arises from 
the finite size of the system and we should not take that 
into account. To fit the data we have used a scaled ver- 
sion of the functional form obtained for the free gas, 

G{a) = " ° (30) 

V2n^J 1 + ^1 + [2aia]2 

where ap and ai are fitting parameters, representing lin- 
ear scaling of the conductance and energy axes. As we see 
from Figure|31the fit is very good over the crucial middle 
range of frequencies. For a > 3 it becomes worse but 
we already know that this difference can be attributed to 
convergence problems with number of K points due to 
the cusp in xi^: ^')- When looking at a — s- limit of our 



7 



fitted expressions (coefficients qq in the figure ISJ we see 
the values bracket the exact value with an relative error 
of roughly 10%. Averaging the result for even and odd 
states gives a smaller and more convergent error (0.5%). 

This calculation is intended to be illustrative, yet it 
shows that the formula H17|l can, in principle, lead to 
the correct answer even when applied to system with 
"zero" boundary conditions. We particularly believe that 
the use of periodic boundary conditions will significantly 
improve the performance of this approach for the main 
problem of the zero boundary condition is clearly the 
one-way transfer of charge from one side to the other, 
which introduces finite-system errors after a relatively 
short time, and therefore spoils G(a) between a = and 
a relatively large value. This gross effect is not present 
for a periodic system. At the same time, periodic bound- 
ary conditions remove the unwanted artificial oscillation 
in G with the number of occupied states, arising from the 
even/odd symmetry as discussed above. 



VI. CONCLUSIONS 

We have developed a unifying point of view of the po- 
larization, the nonlocal conductivity and the conductance 
which supplies a steady-state transport characteristics of 
any system. We have shown that the weight of the Drude 
singularity at zero frequency of the nonlocal conductivity, 
when considered in reciprocal space, directly corresponds 
to the conductance of the system to which we apply a 
field with a nonzero overall drop in potential. We have 
identified a simple relation between conductivity and po- 
larization for the case of a system under unidirectional 
external field that eventually led us to an simple formula 
for conductance, expressed through the polarization of 
the system at small imaginary frequency. Expressed in 
terms of polarization, it turned out to be possible to ad- 
dress the self-consistent field that contributes to the drop 
in total potential used for definition of conductance. We 
have shown that the formula for its evaluation remains 
formally intact but instead of the polarization function 
we need to supply a "transport part" of the polarization. 
The latter is identical to the irreducible polarization in 
ID, but differs from it in general 2D or 3D systems except 
where perfect translational invariance exists perpendicu- 
lar to the flow of current. 

Finally we have demonstrated that our expression for 
conductance in terms of polarization can be used for the 
convenient numerical evaluation of the conductance for 
systems without imposing specific boundary conditions 
in the form of scattering states. This formulation is di- 
rectly suitable for inclusion of many-body or inelastic ef- 
fects, since it is based on the polarization function for 
which approximations that include these complications 
are very well developed. 
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Appendix 

A. Derivation of equation J17I) 

Let g{q,q') be the inverse operator to the right hand 
side in for G = 0. 

J dkg{q, k) (^5{k - q') - ^x{k, ?')) - S{q - q'), (31) 

so that E{q) — J dkg{q,k)E^{k). According to lO we 
have 

Iiq)^jdkj'^^^^gW,k)E\kl (32) 

which motivates us to define 

X\q.q)^q j dk^^g{k^q'). (33) 

Multiplying (|31|l with q'x{q'\ q)/q from left and integrat- 
ing over q we have 

X\q", q') -^ttJ dk^^i^xik, q') = x{q" ,q'), (34) 

comparing term-by-term this Dyson-like equation with 
the Dyson equation for the reducible polarization x{1^ q') 
we arrive at the equation (|18|) which omits the K = 
term from the sum. Substitution of H33|) into (|32|l give 
immediately the result (|17|l which concludes the stated 
results in section Hill 

B. Derivation of equation (I25II using formula II17B 

The polarization of non-interacting electrons is given 
by expression^^ 

X{q,q';^a) = ^(»|e-^-|j)(j|e-'^'"|z) . ^' ~ , (35) 

la ~ Ci + ei 

where |i) are eigenstates of the electronic Hamiltonian 
with eigenenergy e.i, rii is its occupation factor and 

{i\e-^nj)= J dxmx)e-''^^x\j). (36) 

Using formula (|10|l we therefore have 

ij 

(37) 
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The first two matrix elements are, after integrations, 
complex conjugate to each other and therefore their prod- 
uct is real. The last fraction in 1)3 7() is not state- but 
energy-dependent and its real part, which is only needed, 
is given by 



dn{ei 



■{ei 



(38) 



We have already used the fact that in the limit of our 
interest [a 0) the first factor will be sharply peaked 
and therefore we can use linear Taylor expansion of n(ej) 
around e^. However, the expression will be non-zero only 
if the factor (e^ — e^)^ will be compensated by the energy 
dependence of the matrix elements, which we will confirm 
in following paragraphs. 

In the next step we interchange the order of integra- 
tions in (|37(l so that we directly evaluate 



dq 



2-K q- 



(39) 



where we interpret the singularity to be at g = —ie. 
This is in fact arbitrary since, as we have pointed out 
in the section ITlll ^(q, g';ti>) = qq'f{q,q'). An equally 
valid choice of singularity at q = +ie leads to identical 
results. 

We can now finally turn to evaluation of the matrix 
elements. The asymptotic character of the wavefunctions 
has the well known form 



t>k,R{x) 

h.L{x) 



1 f e'''^ + rfec"*''^ a; < 
ifee''^^ X > 

ike"''"'' 



/2tt 
1 



X < 



(40) 
(41) 



The final integrals we need to do have the form 



{k,R/L\k\R/L)= f e'-rk,R/L{x)'Pk',R/L{x)dx, 

J — 00 

(42) 

and of there we need to keep only those that are singular 
~ jTzp" since in H38|l we need to compensate the factor 
(efc — Ck'Y = \(k + k'Y{k — k'y. Foreseeing the delta- 
function character with respect to — Ck' and — ep 
appearing in (|38|l we can directly use tk — t^' — tkp and 
fk = "fk' = fkp I transmission and reflection proba- 
bility amplitudes at the Fermi energy respectively. It is 
now easy to see that the singular terms are 



{k,R\k',R) ^ -~{k,L\k\L) 
{k,R\k',L) = ik,L\k',R)* 



, 1 1^1^ 
2tt k-k' 
1 r*i 
2TTk' -k' 



(43) 
(44) 



Using these forms, together with H38|) and the form zero 
temperature limit that dn{e)/de = 5{e — ep) we directly 
get 



the celebrated Landauer formula. 



\tl 
27r' 



(45) 
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